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Abstract — The entire case study aims to present a procedure for synthesis of discrete steady -state error free modal state 
controller, based on predefined pole placement area and measurable state variables , suggesting a solution for the so called 
pole placement problem (PPP) ( known also as pole assignment problem - PAP). It will introduce a research on how the 
suggested solution will result in modeling of dual-mass DC electromechanical plant. 
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I. Introduction 

One of the most frequently used approaches to dynamic response analysis of linear time -invariant systems modeling, is so 
called modal control, which is related to the assignment of the closed -loop poles to an arbitrary predefined set of points in the 
complex plane. For continuous -time systems the pole placement could be assigned to predefined values or into areas which 
should be distributed in the left half plane of the complex plane. In terms of discrete -time systems areas desired must be 
within the unit circle in the z-plane. Regardless of the method chosen for the pole placement there are two different 
approaches depending on the information used in control processing: output feedback pole placement and state feedback pole 
placement. 

The entire case study will develop the second approach state feedback pole placement for which are known two possible 
solutions - all state space variables are measurable and known and the second option is related to information only partially 
available for the state space variables which is a little bit more complex task in terms of analysis and realization of modal 
control methods. Thus it will be a matter of further investigations. 

II. Formulation of the method used, problem statement and preliminaries 

The modal control will be based on the method related to the second approach of the pole placement procedure - predefined 
pole placement area which in this case in particular will be chosen to be elliptical area of distribution in the unit circle. 

The object taken into consideration for modal control is a dual-mass DC electromechanical plant [2,3]. The block diagram of 
this is presented on Fig.l 



figure 1. Block-diagram of the plant under consideration 

2.1 Continuous state space model composition 

Linear continuous -time system state space model taken into consideration is given as follows: 

x = Ax + bu + v (i) 

y = c T x + du 
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where x is an ^-dimensional state vector, u is an r-dimensional input vector, and y is an m - dimensional output vector, A, b, c, 
d are the system matrices of appropriate sizes and v is the disturbance input matrix. 


Based on the practical point of view a set of state variables is selected, corresponding to real physical parameters, as listed 
below: - armature current i a ; x 2 - rotor speed x% - engine elastic torque M u and x 4 - working/actuating machine 

speed (plant output). According to the block diagram shown on Fig.l the state space variables could be presented by the 
following equations: 


*1 = tV ( k PC U ~ C<PX 2) - 

X 2 =J (cQ>x 1 -X 3 - M lt ) 
x 3 = c 12 (x 2 - x 4 ) 

X4 = — (*3 — M l2 ) 

Therefore matrices A, b,c,d and v of equation (1) are: 
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2.2 Discrete state space model derivation 

Similar to the continuous system the discrete one could be also presented through a model in state space using the form given 
below: 


x{k + 1) = A*x{k ) + b*u(k) + v*(k ) ( 4 ) 

y(fe) = c* T x(k ) + d*u(k) 


For the system, given with equation (4), we can consider the following pole assignment problem (PAP) solution: 
determination of the feedback matrix /“ using the same method for the closed-loop state matrix A“ — has its 
eigenvalues distributed inside the desired area (into z plane). If we assume that is an elliptical area symmetric 
according the real axis, then this will guarantee that the elements of matrix /“ are real numbers. The desired area is given by 
the following equation: 

zOi, r 2 , = n + r ± cos(a)T ) + jr z sin{biT ) = [(^ + r 2 + 2^^ + (?i - r 2 )] (5) 

£S new 

where: r t -major semi-axis; r 2 - minor semi-axis; u - distance from the center to the origin; T ^ - discretization period. 

Immediate response is required for the closed-loop system even if input disturbance occur. Placing the poles in 
implementing appropriate state feedback leads to more or less overdamped time responses of the closed -loop system. The 
result will be settling time lower than the plant and steady-state error free. 

III. Results 

Algorithm for identifying a desired feedback matrix /“ - to satisfy the requirements all the closed-loop poles to be distributed 
on and into the defined area given by equation (5) 

eig(.A x )=eig(Aj') 

Step 1 : Find state matrix of the discrete open-loop system in three diagonal band form A > A j- 
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Step 2: Determination of ‘working’ matrix A Tw through correction of the elements of A ? : A? > A 

Step 3 : Determination of eigenvalues of A^^ t : A* — &ig (yljapt). 

Step 4\ Determination of eigenvalues set on an d into the unit circle: = -—t-., i = Ln. 

Step 5: Ensure the desire eigenvalues to be set on and into the defined area : 

= 7F— [Ol + r 2)KL,l + 2 ^newi + Ol ~ r 2 ~)\ « = 

z A newi 

Step 6 : Determination of desired polynomial _ poly ^ ) ^ ) ' 

In the algorithm shown above A? is three-diagonal band matrix found through the elements of the state matrix A u determined 
in accordance with the requirement for equality of the eigenvalues ( Step 1). Conclusion about system stability could be done 
based on the signs and values of its elements [1,4]. Correction of values of the elements is performed in a manner that 
eigenvalues of A? (eigiAjJ) will be positioned on and in the left side of imaginary axis ( Step 2). There is a number of 
options available for use however the procedure could be reduced only to meet the requirement for positive signs of all the 
sub-diagonal elements of >1^ and also to have one complex conjugated pair poles laying on border. Fulfillment of this 
requirement corresponds to a zero element in the matrix row before the last one. The presence of one real pole on the region 
border corresponds to existence of one zero element on the last matrix row. Desired characteristic polynomial fejtz) 
(determined in Step 6) contains information about the required feedback matrix /“ which could be obtained using the 
following equation: 

fi = <4i ~ <4 = ^n (6) 

where are desired polynomial coefficients and a* are initial polynomial coefficients. 

For the closed-loop system with modal controller the prescribed pole position is ensured initially by input = u - f~x. 
Even if the algorithm given above provides the desired pole placement there is still an essential disadvantage - resulting 
closed-loop system is still not steady state error free according to the output oj 2 . There is a general solution of this problem 
and it could be solved if the input signal is determined by the following dependency: k* X[(jc) — xQi) where 
4- 1) = — c^xO c). Once it is implemented the new state space variable is an output of discrete integrator with 

coefficient k *. The result will be equivalent discrete five state space variables closed-loop system given by the below 
equation: 

l.i A .c .q- n i'-[’ -V 5^12 ii? t k j T G $(jc} T Af jWj Cit 3 

where x = \ x \* = -f ■■ ^ 0 0 0 0] T . ]■■ 

k;- b = \ki -A" -fi -Fi — - matrix of unknown state feedback coefficients and integrator coefficient k*. 

IV. Simulation 

Dual-mass DC electromechanical plant shown on Fig.l is operative through DC motor with the following parameters [2,3]: 

Z/iunn = 220 ^ P nom = 0 = 1000 tr/mm; = 2 A\u - 0 - 10F; 

R a = 20,832811; T a = 0,01265; = l,7021fa; c u — l,7Q2A r m: 

k p€ = 22 ; }\= 0,042 kgm 2 ', J 2 — 0,021 kgm 2 \ M [L = G,3404A r m; Af^ = 3,063 6Nm 

In accordance with equation (3) for continuous -time state space system model we calculate: 
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The calculations are made in MATLAB environment and the discrete model of the system is formulated based on the 
continuous model of equation (1), considering (2) and (3). 


Discrete-time state space model could be obtained using MATLAB function for continuous to discrete model transition - 
c2d(.). Discretization period is defined Tq. = 0,1s which is chosen based on requirement Tq, < (0,05 4* 0,l)t ff . In this case t € 

(forecasted control time) is set to be equal to the expectation for the continuous -time closed-loop system [3]. Considering 
these assumptions for the discrete-time state-space model (4) is calculated: 
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As per the discrete-time state-space model the plant system is fully controllable since the matrix of controllability is 
identified as one of full rank, i.e. rank (Q u y ) = 4. 

The input signal may vary within the range u 1 10 r p J which corresponds to variations of the output voltage of the power 
converter (PC) within the range 0 220 V. The current experiment is performed using an input signal n = 9T for all 

experimental data shown on Fig.2. The disturbance = 3,0636 Nm is applied at the time when the plant processes are 
settled already (at the 8 th second of the experiment start time). 



Ms] t,[s] 


FIGURE 2 . START PROCESSES OF THE PLANT 
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Following the equation (7) five state space variables plant matrices and ^spEcrct and the closed-loop matrices 

G" and are as follows: 
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The parameters of area are assumed as: major semi-axis = 0,3, minor semi-axis r 2 =0,15, distance from the center of 
the elliptic area to the origin & = 0,5 . 

For isolating the elements of the state feedback matrix is used algorithm shown above. The defined algorithm is applied 
to in order to obtain the coefficients a^. Coefficients a* are determined by As a result the state feedback 

matrix is calculated as follows: Iff p = .0,5088 0,1267 0,26 0,4669 —0,0587]. 

The transient processes for the state-space variables are shown on Fig. 3 for the plant and for the closed-loop system: input 
signals u = 9V and g = 917. Disturbance M[ 2 = 3,Q636A r m is applied at the 8 th second. 

For alteration range unification of output variables -T liJ3 [ HlT2 t and ^^coMid-EMp is implemented scale factor with coefficient 

k* = 10 , 21 . 




figure 3. Plant and close-loop state variable start processes 
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The block diagram of the closed loop systems is shown on Fig.4, which is transformed in Simulink/MATLAB modeling 
utilities. The resulted model is used for simulation of the working processes in the system to support the analysis and verify 
the ability of the modal controller to work in combination with the rest system components. 



figure 4. Block diagram of the close-loop system 
V. General conclusions 


• The method of pole assignment is adequately and successfully used for synthesis of discrete model in MATLAB 
environment of modal state controller. 

• The algorithm suggested is fully applicable for modeling and experimental research of the parameters of discrete modal 
steady-state error free state controller with measurable state variables. 

VI. Future developments 


The entire paper concerns only one of the methods possible for pole assignments thus future research will be developed 
involving the case where not all of the state-space variables are measurable and known using observer to create a special 
response (feedback) to compensate (reconstruct) the variables for which we do not have complete information or are hardly 
adaptable by some reason to the algorithm used in the entire research. 
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